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Abstract 



We present results for the process gg — > "fZg within the standard model at hadron colliders. Due 
to large gluon flux at the LHC and the proposed HE-LHC, this process can have significant cross 
section at these machines. Several thousand such events have already been produced at the LHC 
according to the standard model predictions. This process can also be a background to new physics 
signals. We discuss some of the important features of this process related to the structure of its 
amplitude. We also compute the total cross section at the hadron colliders and give some important 
kinematic distributions. A comparison of our results is made with the corresponding NLO calculation 
using MCFM program. We briefly discuss the production of 775 and 77Z via gluon fusion. The 
O | amplitude of these processes are closely related to the gg — > jZg amplitude. 
<D 

1 Introduction 

(N .... , T ^ , 

The quest for new physics signals is continuing at the Large Hadron Collider (LHC). The collider is 

presently running at 8 TeV center-of-mass energy and more than 10 fb _1 of data has been collected 
iy-^ so far. The standard model (SM) is in excellent agreement with the collected data [T|. The recent 
discovery of a fundamental boson in the mass range of 125-126 GeV might confirm the existence 
of the long sought Higgs boson of the SM |2j. With the LHC running in the multi-TeV range, the 
gluon-gluon scattering processes that at the leading-order proceed via quark-loop diagrams can be 
important. Observation of such processes can be a significant test of the SM predictions. Many of such 
processes have been computed in the past. Processes with two or more electro- weak vector bosons in 
the final state have been studied extensively in the literature [3 SI d El [TQl [TH [131 H31 US] - In 
this paper, we study the ^Zg production via gluon fusion in detail. The process gg — > jZg contributes 
to the cross section of the process pp — > jZj + X at the NNLO in a s . The amplitude of this process, 
which occurs at the one-loop level, is finite and gauge invariant. Partial results for this process were 
reported in |11| I15j . We found that its contribution at the LHC, compared to the corresponding 
tree-level contribution, was about 4-5 % for a given set of kinematic cuts applied on the final state 
particles. This percentage would increase at higher center-of mass energy machines such as HE-LHC. 
In this detailed study, we check the decoupling of the top quark at the amplitude-squared and the total 
cross section level. We quantify the contribution of vector and axial- vector parts of the jZg amplitude 
towards the total cross section. We also study the scale dependence and the effect of choosing specific 
parton distribution set. In addition, we consider the processes gg —> 775 and gg —> 77 Z and briefly 
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discuss the results at the LHC. The process gg —> 773 at the LHC has already been discussed in 
jlOL 111] . We include the heavy top-quark in the loop to see if it decouples at the current energies. The 
NLO QCD corrections to the hadronic processes, pp — > jZj + X and pp — > ^yjZ + X are presented in 
[mE]. I n the next section, we discuss the structure of the full amplitude. The method of calculation 
and numerical checks are presented in the section 3. In the section 4, we present our numerical results. 
We conclude in the last section. 




Figure 1: The prototype diagrams for the processes gg — > gjZ. 



2 The Amplitude 

This one-loop process gg — > ^Zg receives contribution from two main classes of quark-loop diagrams 
— pentagon and box type, as shown in Fig. Q. Box class of diagrams are due to the triple gluon 
vertices and can be further divided into three sub-classes. This sub-classification has its own physical 
importance as we will see later. Other diagrams can be obtained by suitable permutation of external 
legs. For each quark flavor, there are 24 pentagon-type and 3 x 6 (= 18) box-type diagrams. However, 
only half of the 42 diagrams are independent. Using Furry's theorem, we can see that axial-vector 
contributions of the box diagrams add up to zero. Therefore, the box diagrams give only vector 
contribution. Pentagon diagrams give both vector and axial-vector contributions. We work with all 
six quark flavors. Except top-quark, we treat all other quarks as massless. The amplitude for our 
process has very nice structure. This is explained in detail in [12j. The full amplitude, after including 
contributions from all six quark flavors, has following form 
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M.v A are color ordered amplitudes for the vector and axial-vector parts of the amplitude. AY, and A 
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quark and i-quark contributions respectively. It is clear from the structure of the amplitude that the 
vector and axial-vector contributions are separately gauge invariant. Therefore, the cross section of 
the process is an incoherent sum of the vector and axial- vector contributions. 

The process gg — > 775 proceeds via pentagon and box type of quark-loop diagrams and the 
amplitude is purely vector type. On the other hand the process gg —> 77 Z proceed via pentagon 
type quark loop diagrams only and its amplitude is purely axial-vector type. The structure of full 
amplitudes for these processes are 
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Our goal is now to compute My and Ma for the light quarks and the top-quark loops. The 
procedure is outlined in the next section. 



3 Calculation and Numerical Checks 

It would be sufficient to focus on the calculation of gg — > jZg amplitude; other two processes can 
be considered as special cases of this process with appropriate changes in color and coupling factors. 
Our calculation is based on the traditional Feynman diagram method. We write down the prototype 
amplitudes for each class of diagrams using the SM Feynman rules. The amplitude of all other 
diagrams are generated by appropriately permuting the external momenta and polarizations in our 
code. Only one prototype pentagon amplitude is sufficient if we keep Z-boson external leg fixed in all 
permutations^] The quark-loop traces are calculated using FORM |18] in n-dimensions for the vector 
contributions. For the axial-vector contribution, the quark-loop trace is calculated in 4-dimensions. 
Interestingly, it turns out that the amplitude for the top-quark in the loop and a massless quark in 
the loop can be written in such a way that they have same analytical form in terms of tensor integrals. 
The dependence on the top-quark mass enters only through the reduction of tensor integrals. This is 
of considerable importance because it reduces the length and evaluation time of the amplitude with 
the top-quark in the loop. In the case of pentagon-type diagrams, the most complicated integral is 
rank-5 tensor integral (E pup,jS ); while for the box-type diagrams, rank-4 tensor integral (D^ upa ) is the 
most complicated one. Five point tensor and scalar integrals are written in terms of box tensor and 
scalar integrals using 4-dimensional Schouten Identity. The box tensor integrals are reduced into the 
standard scalar integrals - Aq, Bq, Cq and Dq that follow from the reduction scheme developed by 
Oldenborgh and Vermaseren [p5] . These reduction procedures were implemented using Mathematica. 
We then obtain FORTRAN routines that compute various tensor integrals in terms of lower-rank 
tensor integrals and scalar integrals. The required scalar integrals are calculated using OneLoop 
library [21]. Because of very large and complicated expression of the amplitude, we calculate the 
amplitude numerically before squaring it. This requires numerical evaluation of the polarization 

1 This is particularly economical when considering axial vector contributions for massive quark case. 
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vectors of gauge bosons. We choose real basis, instead of helicity basis for the polarization vectors to 
calculate the amplitude. This is to reduce the size of compiled program and the time taken in running 
the code. The three-body phase space is generated using RAMBO [22]. 

The amplitude for gg — > ^/Zg is a one-loop amplitude at the leading order (LO) itself. Though 
individual diagrams may be UV and/or IR divergent, the full amplitude should be both UV as well as 
IR finite. IR divergence is relevant to massless quark cases only. All these singularities are encoded 
in various scalar integrals. To make UV and IR finiteness check on our amplitude we have derived 
all the required scalar integrals analytically following t'Hooft and Veltman [23]. We regulate the UV 
divergence of the scalar integrals using dimensional regularization and infrared singularities by giving 
small mass to the quark, the mass regularization. Following are the details of various checks made on 
our amplitude. 

1. UV Finiteness: 

The tadpole and bubble scalar integrals (Aq and -Bo) are sources of UV singularity in any 
one loop amplitude. For the case of massless internal lines, ^4os don't appear in the tensor 
reduction. For both the massive and massless quark contributions, we have verified that the 
amplitude is UV finite. The amplitude of each pentagon diagram has only UV finite tensor 
integrals. Therefore, each pentagon diagram is UV finite by itself as expected from naive power 
counting. The box diagrams individually are not UV finite. Therefore, the cancellation of the 
divergence in the sum of the box diagrams results is an important check. We find that the three 
classes of box diagrams are separately UV finite. 

2. IR Finiteness: 

The diagrams with massless internal quarks have mass singularities. Even with small quark 
mass (like bottom quark) these diagrams may have large logarithms which should cancel for 
the finiteness of the amplitude. There can be log 2 (mg) and log(m 2 ) types of singular terms. 
We have checked explicitly that such terms are absent from the amplitude. We have verified 
that the IR finiteness holds for the three classes of box amplitudes and pentagon amplitude 
separately. 

3. Ward Identities: 

Certain mathematical identities can be obtained by replacing a polarization vector by its 4- 
momentum in any of the pentagon/box amplitudes. This way a pentagon amplitude can be 
written as a difference of two (reduced) box amplitudes and also a box amplitude can be written 
as a difference of two (reduced) triangle amplitudes. In the appendix, we have written these 
identities for a prototype pentagon and a prototype box amplitude. We have verified them 
numerically. These Ward-Identities are important checks on individual diagrams and can also 
be utilized for a systematic study of numerical instabilities in tensor reduction near exceptional 
phase space points. 

4. Gauge Invariance: 

As we have seen, because of the color structure, the vector and axial- vector part of the amplitudes 
do not interfere and are separately gauge invariant. The vector part of the amplitude has gauge 
invariance with respect to three gluons, the photon and the Z-boson. This has been checked by 
replacing the polarization vector of any of these gauge particles by its momentum (e' M (pi) — > ) 
which makes the amplitude vanish. As one would expect the pentagon and three classes of box 
contributions are separately gauge invariant with respect to the 7 and Z— boson. For each gluon, 
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one of the three classes of box amplitudes is separately gauge invariant and further cancellation 
takes place among pentagon and the other two box contributions. The axial-vector part of the 
amplitude is separately gauge invariant with respect to the all the three gluons and the photon. 
We verify that due to explicit breaking of the chiral symmetry in the presence of quark mass, 
the axial part of the amplitude vanishes on replacing Z-boson polarization by its 4-momentum 
only in the m q — > limit. 

5. Decoupling of heavy quarks: 

As a consistency check we have also verified that the vector and axial-vector parts of the am- 
plitude vanish in the large quark mass limit [25] . This feature of the amplitude is very closely 
related to its UV structure. The decoupling holds for each pentagon amplitude and also for 
each class of box amplitudes. 

Apart from making all above checks we have confirmed that the vector part of the *yZg amplitude 
(ignoring coupling factors) reproduces the 775 amplitude in the Mz — > limit within the working 
precision. 

4 Numerical Results 

Based on the procedure outlined above, we can now compute the cross section and examine various 
features of the processes. As we have already mentioned, we compute amplitudes numerically using 
real polarization vectors for gauge bosons. For the processes gg — > ")Zg and gg —> 77Z, there are 
total 48 polarized amplitudes to be computed for each phase space point; while for the gg — > 775 
process, there are 32 such amplitudes. Given the number of diagrams, the number of polarization 
combination and the length of the amplitude, the computation becomes very time consuming. Each 
phase space point evaluation takes about 1.3 seconds on a single machine that we use. We, therefore, 
run the code in parallel environment using AMCI package, a PVM implementation of the VEGAS 
algorithm \27\ I28j . We have used more than 30 cores to run the code in parallel environment. Still 
it takes more than 12 hours to get suitable cross section which includes both massive and massless 
quark contributions. We use ifort compiler on Intel Xeon CPU 3.20GHz machines. 

We first study the importance of the diagrams with the top-quark in the loop. For the pro- 
cess gg — > jZg, in Fig. Q, we plot the ratio of amplitude-squared for five and six quark flavor 
contributions as a function of the top-quark mass for a given phase space point. The partonic center- 
of-mass energy, yfs is fixed at 8Mz- Vector and axial- vector contributions are plotted separately. 
The nit = 4M^ corresponds to the scale above which the top-quarks in the loop can not go on-shell. 
Various slope changes shown in this plot correspond to possibilities of producing one or more final 
state particles via on-shell tt annihilation. In Fig. ([3]), we see that the contribution of the top-quark 
(m t = 175 GeV) to the total cross section is negligible. We also see a knee in the plot at rat = 
This corresponds to the Z-boson production with ti annihilation. The top-quark decouples at around 
100 GeV. We have, therefore, not included its contribution in our results presented below. The run 
time of our code is also reduced by 50%. It is not surprising that the top-quark decouples at such 
a low value for our processes. This is simply because there are four /five quark propagators in each 
box/pentagon diagram leading to a large power of the top-quark mass in the denominator. 

We can divide our numerical results that are presented in this paper into two categories. We 
first discuss theoretical results, related to the structure of the amplitude, keeping Z-boson on- 
shell. Theoretical results presented below include following kinematic cuts: pip > 30 GeV, p^ Z > 
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20 GeV, \rp' ^\ < 2.5, R(^,j) > 0.6. These cuts, together with 72(7,7) > 0-6, are also used to com- 
pute total cross sections for 775 and 77 Z cases. We have chosen factorization and renormalization 
scales as \i = \ip = fj,R = Ep, the transverse energy of the Z-boson. Results are obtained using 
CTEQ6M parton distribution functions [29]. In Fig. Q, we present dependence of total cross section 
on collider centre-of-mass energy to see the effect of large gluon luminosity at higher energies. We 
have already seen that the vector and axial- vector parts of the amplitude are separately gauge invari- 
ant. Their contributions towards total cross section are also included in the figure. The axial-vector 
contribution is only about 10% of the total cross section; this contribution comes from the pentagon 
class of diagrams only. Although the box-contribution to the cross section is not separately gauge 
invariant with respect to the gluons, it is gauge invariant with respect to the 7 and Z-boson. We find 
that more than 70% of the total cross section is due to the box-amplitude only, see Fig. ([5j). The scale 



variation of the total cross section about the central value \i = Ef is shown in Fig. ([6j). On increasing 
the scale by a factor of 2, the cross section decreases by about 25%; it increases by about 40% on 
decreasing the scale by a factor of 2. These large variations are expected because our calculation is 
effectively LO as far as /i dependence is concerned. We see that the cross section falls as we increase 
the scale \x. This is because an increase in the factorization scale increases the cross section due to 
increase in gluon luminosity; but an increase in the renormalization scale decreases the cross section 
because of the decrease in the value of a s (fi). When we increase both the scales at the same time, the 
effect of the change in renormalization scale is stronger. It leads to an overall decrease in the cross 
section with the increase in the scale \i. 

Next, we come to the discussion on our phenomenological results. These results include various 
kinematic distributions related to final state particles. For phenomenological results, we work in the 
narrow width approximation. We allow Z-boson to decay into two leptons in the phase space. In this 
case the kinematic cuts are, 

p j T > 30 GeV, v T > 15 GeV > Pt > 10 GeV, \rp' l > j \ < 2.5, 
i2( 7 ,i),i2(/,7),i2(/,j),i2(/,/) >0.4. 

For convenience we have chosen the scale \i = \ip = /j,r = Mz- In Fig. (|7|), we give cross section 
variation in the range of 8 to 14 TeV center-of-mass energy using both CTEQ611 and CTEQ6M 
parton distribution functions. These numbers do not include the branching ratio of Z — > l + l~ . In 
particular, the cross sections, with CTEQ611 (CTEQ6M) parton distributions, are 65.4 (53.0) fb 
and 202.4 (154.3) fb at 8 TeV and 14 TeV center-of-mass energies respectively. With these cross 
sections, number of gg — > jZg events can be as large as 20000 at 14 TeV LHC, with 100 fb -1 
integrated luminosity . However, to observe these events, one may have to look at Z — > l + l~ decay 
channel; here I can be an electron/muon. So including the branching ratios, one may expect more 
than 1000 gg — > 7-Z(— >■ l + l~)g events. The transverse momentum and rapidity distributions for 
the final state particles are shown in Figs. ([8 13) at 8 TeV center-of-mass energy. We have given 



normalized distributions as they remain same for different choices of parton distributions and/or 
scales. These distributions are characteristically similar at different collider center-of-mass energies 
but at higher energies contribution coming from high px / rapidity region grows while low pt /rapidity 
region contribution goes down. We note that pp is softer as compared to pip. It is because the 
cross section is dominated by the box class of diagrams and in these diagrams, the gluon is emitted 
as a bremsstrahlung radiation. Due to the same reason, i.e., the gluon is emitted more collinearly, 
the rapidity distribution of the gluon is broader as compared to that of the photon. The lepton-pr 
distribution peaks around Mz /2. On the other hand the rapidity distribution of lepton is more central 
compared to the rf distribution. 
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y/S (TeV) 


p% (GeV) 


a LO (pb) 


a NLO (pb) 


< NLO (fb) 


a NNLO /a NLO (%) 


8 


30 


2.202 


2.560 


46.05 (38.25) 


1.8 (1.5) 




50 


1.144 


1.289 


30.49 (25.61) 


2.4 (2.0) 


14 


30 


4.868 


5.490 


158.72 (124.48) 


2.9 (2.3) 


50 


2.608 


2.803 


109.92 (86.61) 


3.9 (3.1) 


35 


30 


14.973 


17.069 


854.09 (606.07) 


5.0 (3.6) 


50 


8.220 


8.898 


607.35 (438.88) 


6.8 (4.9) 



Table 1: Cross sections for the production of pp —> jZj + X at various collider center-of-mass 
energies. We use CTEQ611 PDF set at LO and CTEQ6M PDF set at NLO. The NNLO predictions 
are with CTEQ611(CTEQ6M) parton distribution. The factorization and renormalization scales are 
set to fj,p = fJ-R = Mz- 



We have also compared results of this NNLO calculation with the LO and NLO predictions for 
pp — > j(Z — > i>v)j + X |16j . The LO and NLO results are obtained using parton-level next-to-leading 
order program MCFM0 The comparison is presented after removing the branching ratios in Table 1. 
The table includes results at three different center-of-mass energies and for two values of the pj! 1 . 
We have included the center-of-mass energy of 35 TeV, as it is proposed for the HE-LHC collider. 
The other kinematic cuts are: p> T > 30 GeV, pf iss > 30 GeV, \r]^ j \ < 2.5, R(j,j) > 0.4. This table 
illustrates two facts - 1) the fraction of NNLO events increases with the increase in p^ min , 2) the 
NNLO process becomes more important as we increase the center-of-mass energy. There is an increase 
in the NNLO fraction with an increase in pj mm because, in the NLO events, photon is emitted from 
a quark line; so a larger pl^, " suppresses the NLO contribution more than the NNLO contribution. 
In figure (14) we have compared the normalized p^-distributions at NLO and NNLO, leading to the 
same conclusion. The importance of the NNLO process is more at higher center-of-mass energy simply 
because of the increase in the gluon-gluon luminosity. 

The variation of the total cross section with collider centre-of-mass energy for 775 and 77 Z pro- 
duction via gluon fusion are given in Figs. (15) and (16) respectively. We reconfirm the importance 
of gg —> 775 at the LHC. Like the gg —> ~fZg process, we have checked that the diagrams with the 
top-quark in the loop make negligible contribution. The cross section at the 14 TeV centre-of-mass 
energy is of the order of 2 pb, leading to several thousand events with even 1 fb -1 of integrated 
luminosity. The amplitude for gg — > jjZ is purely axial-vector type and the coupling factors are 
smaller as compared to the axial-vector part of the jZg amplitude. We find that its cross section 
is very small at the LHC, about 0.02 fb. Therefore, this process is unlikely to be observed at the LHC. 



Like other calculations of our types, we have also faced the issue of numerical instability in our 
calculation for certain phase space points. This is a well known issue with the reduction of one-loop 
tensor integrals of higher rank and higher points. This issue is also faced in the evaluation of scalar 
integrals. We face numerical instabilities in the evaluation of pentagon-tensor integrals. This is related 
to the inaccurate evaluation of Gramm determinants in those phase space regions where the linear 
independency of external momenta (modulo 4-momentum conservation) is compromised. The inverse 
Gramm determinants appear in the reduction of tensor integrals. Near exceptional phase space points 
the Gramm determinants become very small and give rise to numerical problems because of loss of 
precision due to large cancellation. This problem can be handled in several ways. One way is to 

2 http:/ /mcfm. fnal.gov/ 
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use higher precision for tensor reduction and scalar integral calculations. This certainly reduces the 
number of exceptional phase space points but the code becomes enormously slow. Another approach 
could be to use special expressions for tensor reduction for such phase space points. This approach 
is some time better. It is important to mention here that none of these two approaches cure the 
problem of numerical instability completely [26J. A more economic and convenient way to proceed 
in this situation is to judiciously ignore the contributions from such phase space points. This one 
can do because we are not doing precision calculations and exceptional phase space points are very 
special and are unlikely to give a significant contribution to the total cross section. We perform a 
gauge-invariance (GI) test on full amplitude for each phase space point. In practice, we put a GI test 
cut of, 5 = 10~ 6 on amplitude-squared. We ignore all those points which fail to pass this test. With 
this cut on 5 and other kinematic cuts, the number of points ignored are below 2%. However, with 
higher px cut and/or less stringent cut on the 6, the number of such points can drop to a level of 
0.01%. We sample about 0.4-0.5 million phase space points to obtain the numerical results. Given 
the volume of phase space, the number of exceptional phase space points is small and it is reasonable 
to assume that the cross section is not dominated by this region of phase space. We find that our 
result depends on this cut very weakly and remains quite stable over the range of 10 -4 — 10 -12 for 
the choice of cut. This also reflects that the exceptional phase space points are few and make small 
contribution. This can be seen in Fig. (17). We have also checked our results for cross section 
calculation by implementing a set of Ward identity tests and its sensitivity on 5 like cut parameter. 
Though the total cross section is quite stable, various distributions (specially the rapidity distributions 
near the edges) are quite sensitive to variation in 5. Edges of distributions define the region where 
exceptional phase space points may lie. Therefore, inaccuracy of the distributions at the edges of 
the phase space is not surprising. We have seen that the exceptional phase space points may defy 
the Ward identity tests sometimes. This might be related to the form of these identities as shown in 
eqns. (8) and (9) Q One can make these 5-like cuts more stringent to get more reliable distributions. 
One can also identify exceptional phase space points at the level of Gramm determinants which may 
be more economical. Phase space points corresponding to large cancellation in Gramm determinants 
can be ignored without putting stringent cuts on 5 and again leading to more reliable distributions. 
A method to implement this criterion is discussed in [26J. 



5 Conclusions 

In this paper, we have presented the results of the cross section calculations for the processes gg — > 
^yZg, gg — > 775 and gg —> 77Z at the hadron colliders. These processes occur at the one-loop level 
via pentagon and box diagrams. We have made a number of checks on our calculation. We have 
verified the finiteness of the amplitude and have checked its gauge invariance with respect to all the 
gauge particles. We have also checked the decoupling of the heavy quark and some of the Ward 
identities which are listed in the Appendix. For the jZg process, the axial-vector contribution does 
not interfere with the vector contribution and it is separately gauge invariant. 

Results are obtained using the parallel running of the code. We have presented the cross sections 
and various kinematic distributions mainly for the gg — > jZj process. We find that the top-quark 
contribution is negligible; and that the vector coupling contributions dominate over the axial vector 
couplings. The typical cross section for the process gg —> jZg is of the order of 100 — 200 fb. 
A comparison with MCFM suggests the contribution of this process to the NLO cross-section of 
PP ~* iZj is about 3 — 4%. One may be able to enhance this percentage further by choosing an 

3 The actual implementation of Ward identities, say for pentagon diagrams, is in the form of, (1 — M i Ap M -2 ) < S. 
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appropriate set of kinematic cuts. The amplitude for gg — > 775 is purely vector type and its typical 
cross section is about 2 pb. Both these processes should be observable at the LHC. However, the 
process gg — > "yyZ is more suppressed due to it being a purely axial-vector pentagon process and it is 
unlikely to be seen at the LHC. 
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6 Appendix: Ward-Identity for Pentagon and Box Amplitudes 

We wish to write down the Ward-Identity for prototype pentagon and box diagrams shown in Fig. 
([!]). We choose replacing Z polarization vector by its 4-momentum to setup pentagon identities. We 
ignore all the irrelevant factors of color and coupling etc. and consider vector part of the amplitude. 
It can easily be extended for axial- vector part of the amplitude. 

M P (pi,P2,P3,P4; ei, e 2 = P2, es, e 4 , e 5 ) = tr I { 5j — ^ — ^ — p 2 j -^7— 

V ?H234 — '' 4 q — nb q f"12 — ,n q r\ ~ < n q r — <'k 

-tr(i 1 I 1 i ■ 1 i -}-tr (i 1 j 1 j 1 , 1 

V 5 ^1234 — m q 4 ^123 m q ^$l~ rn q 1 ftJ V 5 ^1234 — m q 4 ^123 m q ^ $12 ~ m q 1 ^ — m q 

where we have used, p = (^ 12 — m q ) — — m q ). 



=> Mp(pi,P2,P3,PA;ei,e 2 =P2,e 3 ,e4,£5) = Mb(pi,P2 +P3,P4] ei, e 3 , e 4 , e 5 ) 

-M B {p\ +P2,P3,P4;ei,e 3 ,e 4 ,e 5 ). (8) 

Similarly we can write for the prototype box amplitude, 



M B (pi,P2,P3;ei,e 2 =P2,£3,e45) = Mt(puP2 + P3] eu e 3) e 45 ) 

-Mt(pi +P2,P3;ei, e 3j e4 5 ). 



(9) 



Here, £45 denotes the triple gluon vertex attached to the box/triangle diagram. These Ward- 
Identities are shown diagrammatically in Fig. (18). 
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Figure 2: Decoupling of top quark. Vector and 
axial-vector contributions are shown separately. 



Figure 3: Decoupling of top quark in total cross 
section calculation for gg — > ^yZg. 
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Figure 4: Vector and axial-vector contributions of 
total cross section. 




Figure 5: Contribution of box-amplitude towards 
the total cross section. 
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Figure 6: Variation of the total cross section with 
the scale, [i = /j,r = fj,F ■ 



Figure 7: Dependence of total cross section on col- 
lider center-of-mass energy. 
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Figure 8: Transverse momentum distribution of 7 
at 8 TeV centre-of-mass energy. 
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Figure 9: Rapidity distribution of 7 at 8 TeV 
centre-of-mass energy. 
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Figure 10: Transverse momentum distribution of g 
at 8 TeV centre-of-mass energy. 
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Figure 11: Rapidity distribution of g at 8 TeV 
centre-of-mass energy. 
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Figure 12: Transverse momentum distribution of Figure 13: Rapidity distribution of lepton at 8 TeV 
lepton at 8 TeV centre-of-mass energy. centre-of-mass energy. 
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Figure 14: Comparison of normalized pt distribu- 
tions of 7 at NLO and NNLO. NLO distribution is Figure 15: Variation of the total cross section with 
obtained using MCFM. the centre-of-mass energy for gg — > 775. 




8 10 12 
VS (TeV) 



14 16 18 20 



Figure 16: Variation of the total cross section with 
the centre-of-mass energy for gg — > 77Z. 
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Figure 17: Dependence of total cross section for 
775 production on Gl-test cut parameter, S. 
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